NONFACTORIZABLE EFFECTS IN CHARMLESS B DECAYS AND B MESON LIFETIMES 
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Status of nonfactorizable effects in hadronic charmless B decays is reviewed. Implications of new CLEO measurements on B° — * 
7r+7r~ and B — > n' K are discussed. Nonfactorizable effects due to color ocet 4-quark operators are calculated using renormalization 
group improved QCD sum rules. The resultant B-meson lifetime ratio t(B )/t{Hj) agrees with experiment. 
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1 Generalized Factorization 

The nonlcptonic two-body decays of mesons are con- 
ventionally evaluated under the factorization hypothe- 
sis. In the factorization approach, the decay amplitude 
is expressed in terms of factorizable hadronic matrix el- 
ements multiplied by some combinations of Wilson coef- 
ficient functions. To be more specific, the factorization 
hypothesis assumes that the 3-body hadronic matrix el- 
ement (hih 2 \0\M) for the decay M — » h\h 2 is approx- 
imated as th product of two matrix elements (h±\ Ji/jO) 
and (h 2 \J 2 \M) . However, it is known that this approach 
of naive factorization fails to describe the decays pro- 
ceeding through the (class-I) color-suppressed internal 
WK-emission diagrams, though it is at work for decay 
modes dominated by (class-II) external VF-emission di- 
agrams. This implies that it is necessary to take into 
account nonfactorizable contributions to the decay am- 
plitude in order to render the color suppression of internal 
W-emission ineffective. 

Because there is only one single form factor involved 
in the class-I or class-II decay amplitude of B — > PP, PV 
decays, the effects of nonfactorization can be lumped into 
the effective parameters a\ and a 2 : u 



n,2 



-1,2 



+ C 



2,1 



(i) 



where Xi are nonfactorizable terms and receive main 
contributions from color-octet current operators. Since 
\c\/c2\ 1, it is evident from Eq. (1) that even a small 
amount of nonfactorizable contributions will have a sig- 
nificant effect on the color-suppressed class-II amplitude. 
If xi,2 are universal (i.e. process independent) in charm 
or bottom decays, then we have a generalized factoriza- 
tion scheme in which the decay amplitude is expressed 
in terms of factorizable contributions multiplied by the 
universal effective parameters af s 2 . For B — > VV de- 
cays, this new factorization implies that nonfactorizable 
terms contribute in equal weight to all partial wave am- 
plitudes so that al 2 can be defined. It should be stressed 
that, contrary to the naive one, the improved factoriza- 
tion does incorporate nonfactorizable effects in a process 



independent form. Phenomenological analyses of two- 
body decay data of D and B mesons indicate that while 
the generalized factorization hypothesis in general works 
reasonably well, the effective parameters a1 s 2 do show 
some variation from channel to chaniieL especially for 
the weak decays of charmed mesons BoO. An eminent 
feature emerged from the data analysis is that af is 
negative in charm decay, whereas it b^c^mes positive in 
the two-body decays of the B mesonC 
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Phcnomenologically, it is often to treat the number of 
colors iV c as a free parameter to model the nonfactor- 
izable contribution to hadronic matrix elements and its 
value can be extracted from the data of two-body non- 
leptonic decays. Theoretically, this amounts to defining 
an effective number of colors N° s , called l/£ in a, by 
l/iV c eff = (1/N C ) + x- It is clear from (2) that 



Nf(B -► Dtt) = 1.8 - 2.2 fa 2 . 



(3) 



2 Nonfactorizable Effects in Hadronic Charm- 
less B Decays 

What are the_nonfactorizable effects in hadronic charm- 
less B decaysQ? We note that the effective Wilson coeffi- 
cients appear in the factorizable decay amplitudes in the 
combinations a 2i = 
-eff (i = 1,...,5) 
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cii_i and a 2i _i = c 2i _ x 
W°2i (* = 1 , • ■ ■ ^ 5) . As discussed in Sec. 1. nonfactor- 
izable effects in the decay amplitudes of B — > PP, VP 
can be absorbed into the parameters a? . This amounts 
to replacing N c in a| ff by (N° s )i- Explicitly, 
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It is customary to assume in the literature that (N°- ) i 
(N° ) 2 ■ •• ps (7V° ff )io; that is, the nonfactorizable term 
is usually assumed to behave in the same way in tree and 
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penguin decay amplitudes. A closer investigation shows 
that this is not the case. We have argued inB that nonfac- 
torizable effects in the matrix elements of (V— A)(V + A) 
operators are different from that of (V — A)(V — A) op- 
erators. One reason is that the Fierz transformation of 
the (V — A)(V + A) operators 05,6,7,8 is quite different 
from that of (V — A) (V — A) operators Oi,2,3,4 and Og,io- 
Hence, we will advocate that 



Nf(LL) = (Nf) 12M1Q 
Nf(LR) = (Nf) 



(5) 



and that Nf(LR) ^ Nf(LL). In principle, Nf can 
vary from channel to channel, as in the case of charm 
decay. However, in the energetic two-body B decays, 
N% is expected to be process insensitive as supported 
by data 

2.1 Nonfactorizable effects in spectator amplitudes 

To study N* (LL) in spectator amplitudes, we focus 
on the class-Ill decay modes sensitive to the interfer- 
ence between external and internal VF-emission ampli- 
tudes. Good examples are the class-Ill modes: i? ± — > 

etc. Considering B ± — > 



lott , we find that the branching ratio is sensitive to 
1/N£ S and has the lowest value of order 2 x 10 _6 at 



N% s = 00 and then increases with 1/N° S - 
CLEO measurement yields □ 



B(B~- 



= (i-il85 ±0-2) 



x 10" 



The 1997 
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Consequently, 1/N° S > 0.35 is preferred by the dataQ. 
Because this decay is dominated by tree amplitudes, this 
in turn implies that N° S (LL) < 2.9. If the value of 
N° S (LL) is fixed to be 2, the branching ratio fox. posi- 
tive p, which is preferred by the current analysis t9, will 
be of order (0.9 — 1.0) x 10~ 5 , which is very close to 
the central value of the measured one. Unfortunately, 
the significance of i? ± — > um^ is reduced in the recent 
CLEO analysis and only an upper limit is quoted til: 
B(B ± -> tt ± w) < 2.3 x 10~ 5 . Nevertheless, the cen- 
tral value of B(B ± — > to) remains about the same as 
(6). The fact that (LL) ~ 2 is preferred in charmless 
two-body decays of the B meson is consistent with the 
nonfactorizable term extracted from B — > (D, D*)(w, p) 
decays: Nf(B -> Dir) w 2. Since the energy release 
in the energetic two-body decays B — > lot:, B — > Dtt is 
of the same order of magnitude, it is thus expected that 
Nf(LL)\ B ^^2. 

In analogue to the decays B — > D^*\tt, p), the inter- 
ference effect of spectator amplitudes in class-Ill charm- 
less B decay can be tested by measuring the ratios: 
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The ratios Ri are greater (less) than unity when the inter- 
ference is constructive (destructive) . Hence, a measure- 
ment of Ri (in particular i?3)El, which has the advantage 
of being independent of the Wolfcnstcin parameters p 
and 77, will constitute a very useful test on the effective 
number of colors N* (LL). 

During this conference, CLEO has reportetLthe up- 
dated limits on B° — > tt + tt~ and B~ — > tt~tt°: E3 



2?(-B° -> 7r+7r) < 0.84 x 10~ 5 , 

-► tt^tt ) < 1.6 x 10~ 5 . 



(8) 



In particular, the limit on B° — > tt + tt~ is improved by 
a factor of 2. It appears that this decay provides a 
stringent constraint on the form factor F^ 71 . Irrespec- 
tive of the values of N£ s , the predicted branching ratio 
for B° — > 7T + 7r - will easily exceed the current limit if 
F o ""(0) > 0.30. Note that the decay rate of B° -> tt+tt" 
increases slightly with N° S (LL) as it is dominated by 
the tree coefficient a\. For F^^O) = 0.30, we find 
Nf(LL) < 0.20. 

2.2 Nonfactorizable effects in penguin amplitudes 

The penguin amplitude of the class- VI mode B — > <fiK 
is proportional to the QCD penguin coefficients (a% + 
fl 4 + a s) an d hence sensitive to the variation of N° S (LR) 
since 04 is JV^ -stable, but 0,3 and as are -sensitive. 
Neglecting ^-.annihilation and space-like penguin dia- 
grams, we findEl that N^ S (LR) = 2 k-evidently excluded 
from the present CLEO upper limit O 



B(B ± -> tpK*) < 0.5 x 10" 



(9) 



and that l/Nf(LR) < 0.23 or N^(LR) > 4.3 . A simi- 
lar observation was also made inE3. The branching ratio 
of B — > (f>K*, the average of 4>K*~ and 4>K*° mp4es, is 



also measured recently by CLEO with the result E. 



B(B -> (j)K*) = (l.l±oj ± 0.2) x 10 



i-5 



(10) 



We find that the allowed region for N* (LR) is 
4 > N° S (LR) > 1.4. This is in contradiction to the con- 
straint Nf(LR) > 4.3 derived from B ± -> 0if ± . In 
fact, the factorization approach predicts that T(B — ► 
~ T(B — > 4>K) when the VF-annihilation type of 
contributions is neglected. The current CLEO measure- 
ments (j|) and ([n]) are obviously not consistent with the 
prediction based on factorization. One possibility is that 
generalized factorization is not applicable to B — > VV. 
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Therefore, the discrepancy between B(B — ► <j)K) and 
23(i? — > 0A*) will measure the degree of deviation from 
the generalized factorization that has been applied to 
fl — * <f)K* . It is also possible that the absence of 
B — > </>A events is a downward statistical fluctuation. At 
any rate, in order to clarify this issue and to pin down 
the effective number of colors N° S (LR), we urgently need 
measurements of B — > </>A and £? — > 4>K* , especially the 
neutral modes, with sufficient accuracy. 

2.5 B —> iq'K decays 

The published CLEO results on the decay B — > 77'A 

8(5* -> rj'K^ = (6.5±i;S ± 0.9) x 10~ 5 , 

6(5° -> 77' A" ) = (4.7±|:o ± 0.9) x 10- 5 , (II) 

are several times larger than previous theoretical pre- 
dictions y^yilZI in the range of (1 - 2) jSWpv :t was 
pointed out last year by several authors t§ll3c3 that the 
decay rate of B — > rj'K will get enhanced because of the 
small running strange quark mass at the scale mt> and 
sizable ST/ (3) breaking in the decay constants—ig-.and 
/o. Ironically, it was also realized last year thattata the 
above-mentioned enhancement is partially washed out by 
the anomaly effect in the matrix element of pseudoscalar 
densities, an effect overlooked before. Specifically, 

2 

771 / 

(Vl*ft-|0> = -i^ (12) 

where the QCD anomaly effect is manifested by the de- 
cay constant fijf,. Since ~ ^/*, it is obvious that the 
decay rate of -B — > 77'A induced by the (S - P)(S + P) 
penguin interaction is suppressed by the anomaly term 
in (t/|s75s|0). As a consequence, the net enhancement 
is not large. If we treat N° S (LL) to be the same as 
N° S (LR), as assumed in previous studies, we would ob- 
tain typically B(B ± -> 77' A±) = (2 - 3) X 10~ 5 (see the 
dot-dashed curve in Fig. 1). 

What is the role played by the intrinsic charm con- 
tent of the 77' to B — > rj'K ? It has been advocated 
that the new internal ^-emission contribution coming 
from the Cabibbo-allowed process b — > ccs followed by 
a conversion of the cc pair into the 77' via two gluon 
exchanges is potentially important since its mixing an- 
gle V c bV* s is as large as that of the penguin amplitude 
and yet its Wilson coefficient 0,2 is larger than that 
of penguin operators. The decay constant fZ, defined 
by (0|cqv75c|7/) = ifZqp, has been calculated theoret- 
ically Ejo and extracted phenomenologically from the 
data of J/tp — > 7y c 7, J/tp .^.-j/7 and of the 777 and 
7/7 transition form factors Mr 3 !; it lies in the range - 
2.3 MeV < < -18.4 MeV. The sign of is cru- 
cial for the 77' charm content contribution. For a neg- 
ative /5, its contribution to B — ► rj'K is constructive 



for a 2 > 0. Since a 2 depends strongly on N* (LL), we 
see that the cc — > if mechanism contributes construc- 
tively at l/Nf(LL) > 0.28 where a 2 > 0, whereas it 
contributes destructively at 1/N° S (LL) < 0.28 where 0,2 
becomes negative. In order to explain the abnormally 
large branching ratio of B — > 77' A, an enhancement from 
the cc — * 77' mechanism is certainly welcome in order to 
improve the discrepancy between theory and experiment. 
This provides another strong support for N^ S (LL) as 2. 
If Nf(LL) = Nf(LR), then B(B -> r/'K) will be sup- 
pressed at 1/A c eff < 0.28 and enhanced at 1/A c cff > 0.28 
(see the dashed curve in Fig. 1 for /5 = —15 MeV). If 
the preference for A c cff is 1/A c cff < 0.2 (see e.g. 0), then 
it is quite clear that te contribution from the 7/ charm 
content will make the theoretical prediction even worse 
at small 1/A c cff ! On the contrary, if Nf(LL) « 2, the 
cc admixture in the 77' will always lead to constructive in- 
terference irrespective of the value of N° S (LR) (see the 
solid curve in Fig. 1). 
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Figure 1: The branching ratio of B^ 1 — * rj'K ± as a function 
of l/Nf(LR) with Nf(LL) being fixed at the value of 2 
and 77 = 0.34, p = 0.16. The charm content of the rf with 
f°, = —15 MeV contributes to the solid curve but not to the 
dotted curve. The anomaly contribution to (77'|f75s|0) is in- 
cluded. For comparison, predictions for Nf(LL) = Nf(LR) 
as depicted by the dot-dashed curve with /°, — and dashed 
curve with = —15 MeV is also shown. The solid thick 
lines are the preliminary updated CLEO measurements (13) 
with one sigma errors. 

At this conference we learned that a recent CLEO 
reanalysis of B — > 77'A using a data sample 80% larger 
than in previous studies yields the preliminary results 

B{B ± -> 7/'A ± ) = (7.4±5;| ± 1.0) x 10~ 5 , 
B(B d -> 77'A ) = (5.9±i:| ± 0.9) x 10- 5 , (13) 

suggesting that the original measurements (11) were not 
an upward statistical fluctuation. This favors a slightly 
larger /5 in magnitude. For N^ S (LL) = 2 and /5 = —15 
MeV, which is consistent with all the known theoretical 
and phenomenological constraints, we show in Fig. 1 that 
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B(B ± -> rj'K^ at l/Nf(LR) < 0.2 is enhanced consid- 
erably from (2.5-3) x 10~ 5 to (4.6-5.9) x 1CT 5 . In addi- 
tion to the i]' charm content contribution, N° S (LL) m 2 
leads to constructive interference in the spectator ampli- 
tudes of B — ► r/'K and an enhancement in the term pro- 
portional to 2(a 3 — a 5 )xi BK j2 ' + (as + a 4 — a 5 )xi BK ' v \ 
In the B s system, we Ej find that B s — > 7777', the 
analogue of £?° — * rfK°, has also a large branching ratio 
of order 2 x 1CT 5 . 

i?.^ Summary 

Tree-dominated rare B decays, i? ± — * utt 111 and B° — ► 
7T+7T-, favor a small A c cff (LA), namely Nf(LL) » 2. 
The constraints on N° (LR) derived from the penguin- 
dominated decays -B 1 * 1 — * i^A 1 * 1 and B — > 0A*, which 
tend to be larger than N° (LL), are not consistent with 
each other. Our analysis of B — > 77' A clearly indicates 
that Nf(LL) w 2 is favored and Nf(LR) is preferred 
to be larger. The preliminary updated CLEO measure- 
ments of _B — * 77' X seem to imply that the contribution 
from the 77' charm content is important and serious. 

3 Final-state interactions and B — > w/Y 

The CLEO observationE] of a large branching ratio for 
B ± -► wA ± 

5(5* -> wA*) = (l.5±8:e ± 0.2) x 10~ 5 , (14) 

is difficult to explain at first sight. Its factorizable am- 
plitude is of the form 

A(B~ -> wK~) cx (a 4 + Ra 6 )X^ BuJ ^ K) (15) 
+ (2a 3 + 2a 5 + ia 9 )X( SK ^ + ■ • ■ , 

with i? = — 2rn 2 K /(rnt > rn s ), where ellipses represent for 
contributions from IF-annihilation and space-like pen- 
guin diagrams. It is instructive to compare this decay 
mode closely with B~ — > pK~ 

A{B- -» p°K-) cx (a 4 + i?a 6 )^ (Bp ' K) + • •• • (16) 

Due to the destructive interference between a 4 and a§ 
penguin terms, the branching ratio of B^ — > p°K^ is es- 
timated to be of order 5 x 10 -7 . The question is then why 
is the observed rate of the wK mode much larger than the 
pK mode ? By comparing ( |l5|) with ( |l6|) , it is natural to 
contemplate that the penguin contribution proportional 
to (2d3 + 2d5 + iag) accounts for the large enhancement 
of B^ — > ojK^. However, this is not the case: The coef- 
ficients as and as , whose magnitudes are smaller than a 4 
and a6, are not large enough to accommodate the data 
unless Nf(LR) < 1.1 or Nf(LR) > 20 (see Fig. 9 ofH). 



So far we have neglected three effects in the consid- 
eration of B^ — > uiK^: IF-annihilation, space- like pen- 
guin diagrams and final-state interactions (FSI) . It turns 
out that FSI may play a dominant role for B^ — > uiK^. 
The weak decays B~ — > A"*~7r° via the penguin process 
b — > smt and — > A*°7T _ via 6 — * scZd followed by the 
quark rescattering reactions {A*~7r°, K*°tt~} — > 
contribute constructively to — > wif" (see Fig. 2), 
but destructively to £>~ — > pK~ . Since the branching 
ratios for — > K*~ir° and K*°ir~ are large, of order 
(0.5 — 0.8) x 10 -5 , it is conceivable that a large branching 
ratio for B^ — > can be achieved from FSI via inelas- 
tic scattering. Moreover, if FSI dominate, it is expected 
that B(B ± -> wK*) « (1 + V2) 2 B{B° -> uK°). 




Figure 2: Contributions to B~ — > K~ui from final-state inter- 
actions via the weak decays _B — — > K*~ir° and -B - — > K*°-k~ 
followed by quark rescattering. 



4 Nonspectator Effects and B Meson Lifetimes 

In the heavy quark limit, all bottom hadrons have the 
same lifetimes, a well-known result in the parton pic- 
ture. With the advent of heavy quark effective theory 
and the OPE approach for the analysis of inclusive weak 
decays, it is realized that the first nonperturbative cor- 
rection to bottom hadron lifetimes starts at order 1/m 2 
and it is model independent. However, the 1/m 2 correc- 
tions are small and essentially canceled out in the lifetime 
ratios. The nonspectator effects such as IF-exchange 
and Pauli interference due to four-quark interactions are 
of order I/tuq, but their contributions can be poten- 
tially significant due to a phase-space enhancement by 
a factor of 16ir 2 . As a result, the lifetime differences of 
heavy hadrons come mainly from the above-mentioned 
nonspectator effects. 

The four-quark operators relevant to inclusive non- 
leptonic B decays are 

Ov-A = b L ^^q L qhl^bh , 



4 



S -p = b R q L q L b R , 

T v -A = bLlufqL qL^t a b L , 

T S - P = b R t a q L q L t a b R . (17) 

From which one can followEI to define four hadronic pa- 
rameters Bi,B2,£x, £2 relevant to our purposes: 



1 



2m 1 



-(B\0 v - A( s- P) \B) = 



f B m B 



_ JB 



B\(2 



1 ™, 



2m f 



(B\T V -a{s-p)\B 



■7T\ _ /l m B 



-£l(2) 



(18) 



Under the factorization approximation, Bi = 1 and £$ = 
0. To the order of 1/mL the i?-hadron lifetime ratios are 
given by 



t(B- 



r(B° g ) 
r(B° d ) 



1 + 0.043^1 + 0.0006S 2 - 0.6l£i + Q.17e : 



= 1 



+ (-1.7 x 10~ 5 B t + 1.9 x 10~ 5 B 2 
0.0044£i + 0.0050 £ 2 ). 



(19) 



It is clear that even a small deviation from the factoriza- 
tion approximation e% = can have a sizable impact on 
the lifetime ratios. 

We have derived in heavy quark effective theorx the 
renormalization-group improved QCD sum rules E£l for 
theJiadronic parameters B\, S 2) £i> and £ 2 - The results 
areEfl 

Bx{m b ) = 0.96 ± 0.04, B 2 (m b ) = 0.95 ± 0.02, 
ei(m b ) = -0.14 ±0.01, £ 2 (m b ) = -0.08 ± 0.01,(20) 

to the zeroth order in l/m&. The resultant B-meson 
lifetime ratios are T{B~)/r{Bd) = 1.11 ± 0.02 and 
t(BAIt(B,i) ~ 1, to be compared with the world aver- 
agesll: r(S-)/r(S d ) = 1.07 ±0.03 and t(B s )/t(B c1 ) = 
0.94 ± 0.04. Therefore, our prediction for r(B~)/r(B d ) 
agrees with experiment. 
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